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Abstract 

Let P be a set of 2n points in the plane, and let Mc (resp., -Mnc) denote a bottleneck 
matching (resp., a bottleneck non-crossing matching) of P. We study the problem of 
computing Mnc- We first prove that the problem is NP-hard and does not admit 
a PTAS. Then, we present an 0(n 15 log ' 5 n)-time algorithm that computes a non- 
crossing matching M of P, such that bn(M) < 2VT0 • 6n(M N c), where bn{M) is the 
length of a longest edge in M. An interesting implication of our construction is that 
bn(M NC )/bn{M c ) < 2\/l0. 



*Work by A.K. Abu-Affash was partially supported by a fellowship for doctoral students from the 
Planning & Budgeting Committee of the Israel Council for Higher Education, and by a scholarship for 
advanced studies from the Israel Ministry of Science and Technology. Work by A.K. Abu-AfTash and Y. 
Trabelsi was partially supported by the Lynn and William Frankel Center for Computer Sciences. Work 
by P. Carmi was partially supported by grant 2240-2100.6/2009 from the German Israeli Foundation for 
scientific research and development, and by grant 87212211 from the Israel Science Foundation. Work by 
M. Katz was partially supported by grant 1045/10 from the Israel Science Foundation. Work by M. Katz 
and P. Carmi was partially supported by grant 2010074 from the United States - Israel Binational Science 
Foundation. 



1 Introduction 



Let P be a set of 2n points in the plane. A perfect matching M of P is a perfect matching 
in the complete Euclidean graph induced by P. Let bn(M) denote the length of a longest 
edge of M. A bottleneck matching Mq of P is a perfect matching of P that minimizes bn(-). 
A non-crossing matching of P is a perfect matching whose edges are pairwise disjoint. In 
this paper, we study the problem of computing a bottleneck non-crossing matching of P; 
that is, a non-crossing matching Mnc of P that minimizes bn(-), where only non-crossing 
matchings of P are being considered. 

The non-crossing requirement is quite natural, and indeed many researches have con- 
sidered geometric problems dealing with crossing-free configurations in the plane; see, 
e.g. p|-[5||14|. In particular, (bottleneck) non-crossing matching is especially important in 
the context of layout of VLSI circuits [II] and operations research. It is easy to see that 
there always exists a non-crossing matching of P (e.g., match each point with the first point 
to its right). Actually, any minimum weight matching of P is non-crossing. However, as 
shown in Figure [TJ which is borrowed from [6] , the length of a longest edge of a minimum 
weight matching can be much larger than that of a bottleneck non-crossing matching. 



(a) (b) (c) 

Figure 1: (a) The set P. (b) A bottleneck non-crossing matching M\ of P. (c) A minimum 
weight matching M2 of P. Notice that bn^M^) /bn(M\) — > n — 2. 



1.1 Related work 



Matching problems play an important role in graph theory, and thus have been studied ex- 
tensively, see [l3] . The various matching algorithms developed for general weighted graphs 
of course apply in our setting. However, it turns out that one can do better in the case 
of points in the plane. Vaidya 15 presented an 0(ra 5 / 2 log 4 ra)-time algorithm for com- 
puting a minimum weight matching, based on Edmonds' 0(n 3 ) algorithm. Subsequently, 
Varadarajan | l6j described an 0(n 3 / 2 log 5 n)-time algorithm for this problem. For the 
bipartite version, Vaidya 15 presented an 0(n 5//2 log n)-time algorithm and Agarwal et 
al. [I] presented an 0(ra 2+£ )-time algorithm; both algorithms are based on the Hungarian 
method 13 . As for bottleneck matching, Chang et al. 17] obtained an 0(ra 3//2 log 
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time algorithm for computing a bottleneck matching, by proving that such a matching is 
contained in the 17RNG. Efrat and Katz [9] extended this result to higher dimensions. For 
the bipartite version, Efrat et al. [8] presented an 0(n 3 / 2 logn)-time algorithm. Algorithms 
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for other kinds of matchings, as well as approximation algorithms for the problems above, 
have also been developed. 

Self-crossing configurations are often undesirable and might even imply an error con- 
dition; for example, a potential collision between moving objects, or inconsistency in the 
layout of a circuit. Many of the structures studied in computational geometry are non- 
crossing, for instance, minimum spanning tree, minimum weight matching, Voronoi dia- 
gram, etc. Jansen and Woeginger [10] proved that deciding whether there exists a non- 
crossing matching of a set of points with integer coordinates, such that all edges are of 
length exactly d, for a given integer d > 2, is NP-complete. Carlsson and Armbruster [6] 
proved that the bipartite version of the bottleneck non-crossing matching problem is NP- 
hard. Alon et al. [4] considered the problem of computing the longest (i.e., maximum 
weight) non-crossing matching of a set of points in the plane. They presented an approx- 
imation algorithm that computes a non-crossing matching of length at least 2/ir of the 
length of the longest non-crossing matching. Aloupis et al. [5] considered the problem of 
finding a non-crossing matching between points and geometric objects in the plane. See 
also [2j[3|[l4] for results related to non-crossing matching. 

1.2 Our results 

We begin by proving (in Section 2) that the problem of computing Mnc is NP-hard. 
Our proof is based on a reduction from the planar 3-SAT problem, and is influenced 
by the proof of Carlsson and Armbruster mentioned above. As a corollary we obtain 
that the problem does not admit a PTAS. Next, in Section 3, we present an algorithm 
for converting any (crossing) matching M x into a non-crossing matching M = , such that 
bn(M = ) < 2\/l0 • bn(M x ). The algorithm consists of two stages: converting M x into an 
intermediate (crossing) matching M' x with some desirable properties, and using M' x as a 
"template" for the construction of M = . The algorithm implies that (i) Mmc/Mc < 2\/l0, 
and (ii) one can compute, in 0(n 3 / 2 log 1 / 2 n)-time, a non-crossing matching M, such that 
bn(M) < 2^/W-bn(M^c)- We are not aware of any previous constant-factor approximation 
algorithm for the problem of computing Mnc- I n the full version of this paper, we also 
present an 0(n 3 )-algorithm, based on dynamic programming, for computing Mnc when 
the points of P are in convex position. 

2 Hardness Proof 

In this section, we prove the following theorem. Our proof is influenced by the proof of 
Carlsson and Armbruster for the bipartite version [6j. 

Theorem 2.1. Let P be a set of In points in the plane. Then, computing a bottleneck 
non-crossing matching of P is NP-hard. 
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Proof: The proof is based on a reduction from the planar 3-SAT problem. Given a 3-CNF 

formula F with n variables X = {x\, X2, ■ ■ ■ , x n } and m clauses Y = {Ci, C2, ■ ■ • , C m }, let 

Gp = (V, E) be the graph of F, i.e., V = XUY and E = {(a;,, Cj) : Xj appears in Cj (either negated or unnegat 

If Gf is planar, then F is called a planar 3-CNF formula. The planar 3-SAT problem is 

to determine whether a given planar 3-CNF formula F is satisfiable; the problem is NP- 

complete fl2] . 

Let F be a planar 3-SAT formula. We construct, in polynomial time, a set P of points 
in the plane, such that F is satisfiable if and only if there exists a non-crossing matching of 
P with bottleneck 1. Consider the graph Gp. It is well known that Gf can be embedded 
in the plane in polynomial time. 
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Figure 2: The true and the false matchings of a variable x. 

Variables. Each variable x% G X is mapped to a circuit Sj of an even number of 
points, where the distance between any two adjacent points is 1. We mark the points of Si 
alternately by p a and p e . Each circuit Sj can be partitioned into pairs of adjacent points 
in two ways. We arbitrarily associate one of them with the assignment x% = T and call 
it the "true matching", and the other with the assignment x% = F and call it the "false 
matching"; see Figure [2] Thus, the value of X{ will determine the matching on Sj, and vice 
versa. 
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Figure 3: The hexagon hj corresponding to the clause Cj. 

Clauses. Each clause Cj £ Y is mapped to a hexagonal component hj of points, where 
the distance between any two adjacent points is 1, as shown in Figure [3j 

Edges. Each edge (xj, Cj) between a variable Xi and a clause Cj is mapped to a path 
lij of an even number of points, where the distance between any two adjacent points is 
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1, that begins and ends at two different points of Si and intersects the hexagon hj, as 
described below. 



C k = (xl V ■ V ■) Cj = (xi V ■ V ■) 
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Figure 4: The paths between Xi and the clauses Cj and Cfc. Since Xj appears unnegated in 
Cj and negated in Ck, h,j begins at p e and ends at p a , and begins at p and ends at p e . 

Assume that Xi appears unnegated in clause Cj and negated in clause C\~, i.e., Cj = 
(xi V ■ V ■) and Cj. = (a^ V • V •). Consider the points of Sj in clockwise order. For the 
unnegated instance, we select 4 consecutive points p e ,Po,Pe,Po of Sj, connect to the first 
and last of these points, and remove the middle two points; see Figure Qa). And, for the 
negated instance, we select 4 consecutive points Po,Pe,Po,Pe, connect & to the first and 
last of these points, and remove the middle two points. Notice that, as for circuits, a path 
can be partitioned into pairs of adjacent points in two ways. The value of Xi and whether 
or not it appears negated in the clause will determine which of these two matchings will 
be the matching on the path, see below. 
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Figure 5: The intersection between Cj and the variables appearing in it. 

The intersection between lij and the hexagon hj representing Cj is shown in Figure [5] 
kj intersects hj such that six of its points lie inside hj and the other lie outside hj. We 
would like to prevent situations where a point of lij is matched to a point of hj. To this 
end, we add 3 pairs of points around each of the two junctions involving kj and hj, see 
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squares A and B in Figure [5j This forces the points of kj to be matched to each other. 

Notice that in the true (resp., false) matching of Sj, if Xi appears unnegated (resp., 
negated) in Cj, then the six points of kj in hj are matched to each other (see Figure [4](b)) , 
and, if x% appears negated (resp., unnegated) in Cj, then the two extreme points among 
these six points are matched to points outside hj. 
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Figure 6: Cj = (x, V V x;). (a) Assigning values to the variables Xi,Xk,xi so that the 
value of each of the corresponding literals is false, ensures that a matching on hj with 
bottleneck 1 does not exist, (b) Assigning value to the variables so that the value of at 
least one of the literals is true, ensures that such a matching does exist. 

To see the correctness of the reduction, consider Figure [HJ First observe that there 
exists a matching on hj with bottleneck 1 if, and only if, for at least one of the three paths 
intersecting hj, its six points lying in hj are matched to each other (i.e., the value of the 
corresponding literal is true). Indeed, in Figure [6]^a) , the value of each of the three literals 
of Cj is false, ensuring that a matching on hj with bottleneck 1 does not exist. And, in 
Figure [6^b), the value of one of the literals (a?*) of Cj is true, ensuring that such a matching 
does exist. Notice that if the value of x~k and/or x\ were also true, then one could still use 
the same matching on hj. 

Conversely, assume there exists a matching with bottleneck 1 and consider the truth 
assignment implied by the matchings on the circuits. We need to verify that each of 
the clauses is satisfied by this truth assignment. Let hj be the hexagon representing Cj. 
Since there exists a matching on hj with bottleneck 1, then, for at least one of the paths 
intersecting hj, its six internal points are necessarily matched to each other. But this 
implies, as noted above, that the value of the corresponding literal, and therefore of Cj, is 
true. 

■ 

Finally, we observe that if one can find a non-crossing matching of bottleneck less than 
3^/2/4, then one can solve the planar 3-SAT problem in polynomial time. This bound is 
obtained from examining the additional points that are added to each hexagon (see squares 
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A and B in Figure [5]). 

Corollary 2.2. The problem of computing a bottleneck non-crossing matching does not 
admit a PTAS, unless P=NP. 

3 Approximation Algorithm 

Let P be a set of 2n points in general position in the plane. The bottleneck of a perfect 
matching M of P, denoted bn(M), is the length of a longest edge of M. Let M x be a 
perfect matching of P. In this section we show how to convert M x into a non-crossing 
perfect matching M= of P, such that bn(M=) < 2 V / 10 • bn(M x ). 

Set <5 = bn(M x ). We begin by laying a grid of edge length 2^/25. W.l.o.g. assume that 
each of the points in P lies in the interior of some grid cell. Consider an edge e of M x . 
Since e is of length at most 5, it is either contained in a single grid cell, or its endpoints lie 
in two adjacent cells (i.e., in two cells sharing a side or only a corner). In the former case, 
we say that e is internal, and in the latter case, we say that e is external. We distinguish 
between two types of external edges: a straight external edge (or s-edge for short) connects 
between a pair of points in two cells that share a side, while a diagonal external edge (or 
d-edge for short) connects between a pair of points in two cells that share only a corner. 
Finally, the degree of a grid cell C, denoted deg(C), is the number of external edges with 
an endpoint in C. 

Our algorithm consists of two stages. In the first stage we convert M x into another 
perfect matching M' x , such that (i) each edge of M' x is either contained in a single cell, or 
connects between a pair of points in two adjacent cells, (ii) for each grid cell C, deg(C) < 4 
and these deg(C) edges connect C to deg(C) of its adjacent cells, and (hi) some additional 
properties hold (see below). In the second stage, we construct the matching M= according 
to M' x , such that there is a one-to-one correspondence between the external edges of M' x 
and the external edges of M=. That is, there exists an edge in M' x connecting between 
two adjacent cells C\ and C2 if and only if there exists such an edge in M=. However, the 
endpoints of an external edge of M= might be different than those of the corresponding 
edge of M' x . 

The second stage itself consists of two parts. In the first part, we consider each non- 
empty grid cell separately. When considering such a cell C, we first determine the at 
most four points that will serve as endpoints of the external edges of C (as dictated by 
M' x ). Next, we construct a non-crossing matching for the remaining points in C. In the 
second part of this stage, we add the external edges between the points that were chosen 
as endpoints for these edges in the first part. 

We now describe each of the stages in detail. 
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3.1 Stage 1 

In this stage we convert M x into M' x . We do it by applying a sequence of reduction rules 
to the current matching, starting with M x . Each of the rules is applied multiple times, 
as long as there is an instance in the current matching to which it can be applied. When 
there are no more such instances, we move to the next rule in the sequence. 



Figure 7: Rule I (left) and Rule II (right). The solid edges are replaced by the dashed 
edges. 

We associate a d-edge connecting between two cells with the corner shared by these 
cells. 

Rule I is applied to a pair of d-edges associated with the same corner and connecting 
between the same pair of cells, see Figure [7} The d-edges are replaced by a pair of internal 
edges. 

Rule II is applied to a pair of d-edges associated with the same corner and connecting 
between different pairs of cells, see Figure [7} The d-edges are replaced by a pair of s-edges. 
Notice that when we are done with Rule I, each corner has at most two d-edges associated 
with it, and if it has two, then they connect between different pairs of cells. Moreover, 
when we are done with Rule II, each corner has at most one d-edge associated with it. 
Finally, since the length of a d-edge is at most <5, any edge created by Rule I or Rule II is 
contained in the disk D${a) of radius 5 centered at the appropriate corner o. 

A d-edge associated with a corner a defines a triangular danger zone in each of the two 
other cells sharing a, see Figure [8^a). A danger zone is semi-open; it does not include the 
hypotenuse. Let S$(a) denote the square that two of its sides are the hypotenuses of the 
danger zones defined by the d-edge associated with o. Notice that if p is a point in cell C 
that is not in the danger zone in C, then one cannot draw a d-edge between C\ and C2, 
with endpoints in the interior of Ss{a), that passes through p. 

Rule III is applied to a d-edge e\ and an edge e2 with an endpoint in a danger zone defined 
by ei, see Figure |8jb-d) . We distinguish between two cases. If ei is a s-edge, then, by the 
claim below, its other endpoint is in one of the cells C\ or C<i- In this case, we replace e\ 
and e2 with an internal edge and a s-edge, see Figure [8^b). If e2 is an internal edge, then 
consider the other endpoint q of e%- If q is not in a danger zone in C defined by another 
d-edge, then replace e\ and e2 with two s-edges, see Figure fflc). If, however, q is in a 
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Figure 8: (a) The danger zones denned by a d-edge. Rule III: (b) &2 is a s-edge. (c) e2 is 
an internal edge and q is not in another danger zone, (d) e2 is an internal edge and q is in 
another danger zone. 



danger zone in C defined by another d-edge e^, then, by the claim below, is associated 
with one of the two corners of C adjacent to the corner a (to which e\ is associated), and 
therefore, either C\ or C2 contains an endpoint of both e\ and e%. Replace e\ and ei with 
two s-edges, such that one of them connects q to the endpoint of e\ that is in the cell that 
also contains an endpoint of es, see Figure [8^d). 

Claim 3.1. Consider Figure^ In an application of Rule III, (i) if e2 is a s-edge, then its 
other endpoint q is in one of the cells C\ or C2, and (ii) if q is in a danger zone defined 
by another d-edge e^ associated with corner b, then ab is a side of C and the danger zone 
containing q is on the side of b containing a. 

Proof. Statement (i) is surely true just before the first application of Rule III, since the 
length of e2 then is at most 5. (Notice that a s-edge created by Rule II cannot have an 
endpoint in a danger zone.) For the same reason, Statement (ii) is surely true just before 
the first application of Rule III. (Notice that if e2 is an internal edge created by Rule I, 
then it is contained in one of the two danger zones defined by e%.) It remains to verify 
that if e2 was created by a previous application of Rule III, then both statements are still 
true. Indeed, if e2 is a s-edge created by a previous application of Rule III, then it had to 
be an application of the type depicted in Figure [8^d) , and the replacement instructions for 
this type ensure that Statement (i) is true. As for Statement (ii), if e2 was created by an 
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application of Rule III, then, since q was an endpoint of a d-edge that was removed by the 
application of Rule III, it is not in a danger zone. ■ 




Figure 9: Rule IV (left) and Rule V (right). 

Let C(p) denote the cell containing point p. 
Rule IV is applied to a d-edge e\ = (pi,qi) and to a s-edge e2 = (j)2,Q2), such that 
C{pi) = C{p2) and C(q\) and C(q2) share a side, see Figure [9} e\ and e2 are replaced by 
an internal edge and a s-edge. 

Rule V is applied to a pair of s-edges e\ = (pi,qi) and &2 = (P2,Q2), such that C{jp\) = 
C(p2) and C(qi) = C(q>z), see Figure[9j e\ and e2 are replaced by a pair of internal edges. 

Let M x be the matching that is obtained after applying Rules I-V. The following lemma 
summarizes some of the properties of M' x ; its proof follows immediately from the discussion 
above. 

Lemma 3.2. M' x has the following properties: 

1. Each edge is either contained in a single cell, or connects between a pair of points in 
two adjacent cells. 

2. A corner has at most one d-edge associated with it. 

3. A d-edge is of length at most 6. 

4- The two danger zones defined by a d-edge e are empty of points of P. 

5. For each grid cell C, deg(C) < 4 and these deg(C) edges connect C to deg(C) of its 
adjacent cells. 

6. If e is a d-edge in M' x connecting between cells C\ and C2, and C is a cell sharing a 
side with both C\ and C2, then there is no s-edge in M' x connecting between C and 
either C\ or C2. 

3.2 Stage 2 

In this stage we construct M = according to M' x . This stage consists of two parts. 
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3.2.1 Part 1: Considering each cell separately 

In this part we consider each non-empty grid cell separately. Let C be a non-empty grid 
cell and set Pc = P n C. We have to determine which of the points in C will serve as 
endpoints of external edges. The rest of the points will serve as endpoints of internal edges. 
We have to consider both types of external edges, d-edges and s-edges. We first consider 
the d-edges, then the s-edges, and, finally, after fixing the endpoints of the external edges, 
we form the internal edges. 

For each corner a of C that has a d-edge (with endpoint in C) associated with it in M' x , 
consider the line through a supporting C and parallel to its appropriate diagonal, and pick 
the point p a in C that is closest to this line as the endpoint (in C) of the corresponding 
d-edge in M = . By Lemma [33] (Property 3), Pc H Ds(a) ^ 0, and therefore the distance 
between p a and a is less than y/25. Moreover, since the length of a side of C is 2y25, each 
of the relevant corners is assigned a point of its own. Observe also that a d-edge in M = 
with endpoint in C will not cross any of the s-edges in M = with endpoint in C. This follows 



from Lemma 3.2 (Property 6). We thus may ignore the points in C that were chosen as 



endpoints of d-edges, and proceed to choose the endpoints of the s-edges. 
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Figure 10: (a) A(p, "up"). Procedure 1: (b) D = {"left", "up", "right"} and \Q\ = 3. (c) 
p A(p', "right") and p is assigned to "up", (d) p £ A(p', "right") and p is assigned to 
"right" ; p" will be assigned to "up" . 

Let A(p,dir) denote the triangle whose corners are p and the two corners of C in 
direction dir, where dir £ {"up", "down", "left", "right"}; see Figure 10 a). If there is only 
one s-edge in M' x with endpoint in C, and its direction is dir, then we pick the extreme 
point in direction dir as the endpoint (in C) of the corresponding s-edge in M=. Assume 
now that there are k, 2 < k < 4, s-edges in M' x with endpoints in C. We pick the k 
endpoints (in C) of the corresponding s-edges in M= according to the recursive procedure 
below. 



Lemma 3.3. The s-edges in M= with endpoints in C do not cross each other. 

Proof. By induction on k. If k = 1, then there is nothing to prove. Assume k > 2 and 
consider Procedure 1. If \Q\ = k, then, for each pi £ Q, the triangle A(pi,diri) is empty, 



see Figure 10 ]b). Therefore, any two s-edges, one with endpoint pi and direction diri and 
another with endpoint pj and direction dirj, i ^ j, do not cross each other. If \Q\ < k, 
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Procedure 1 PickEndpoints (a cell C, a set of directions D = {diri, . . . , dir^}) 
1: Q<-0 

2: for each direction diri 6 D do 

3: let j>j be the extreme point in direction dir^ 

4: Q <r- Q U {p,} 

5: if \Q\ = k then 

6: assign pj to direction dirj, i = 1, . . . , k 
7: Return(Q) 

8: let p £ Q be a point that is the extreme point in two directions dirt, dirj £ D 

9: Q' <- PickEndpoints(C \ {p}, D \ {din}) 
10: let p' G Q' be the point assigned to direction dirj 
11: if p ^ A(p',dirj) then 
12: assign p to direction diri 
13: else 

14: Q' <- PickEndpoints(C \ {p}, D \ {dirj}) 
15: assign p to direction dirj 
16: Return(Q / U {p}) 



then, since k > 2, the directions dirj and dir^ are not a pair of opposite directions. Now, 
if p, the extreme point in directions diri and dirj, is not in A(p',dirj) (Figure 10 c)), 



then \Q' U {p}| = k and each of the corresponding k triangles is empty. If, however, 
p G A(p',dirj) (Figure [To|(d) ) , then let p" be the point assigned to direction dir, L by the 
call to PickEndpoints in Line 14. We claim that p A(p", diri). Indeed, if p were in 
A(p", dir^, then either p" £ A(p', dirj) or p' £ A(p", dir^), contradicting in both cases the 
induction hypothesis for k — 1. ■ 

We are now ready to form the internal edges. Let Pq C be the set of points that 
were chosen as endpoints of external edges, and set Pq = Pc \ Pq ■ We show below that 
P c is contained in the interior of a convex region R C C, such that any external edge with 
endpoint in C does not intersect the interior of R. Hence, if we form the internal edges 
by visiting the points in P c from left to right and matching each odd point with the next 
point in the sequence, then the resulting edges do not cross each other and do not cross 
any of the external edges. 

It remains to define the convex region R. For each endpoint p% of a s-edge, draw a 
line li through pi that is parallel to the side of C crossed by the s-edge. Let hi be the 
half-plane defined by U and not containing the s-edge, and set Ri = hid C. Similarly, for 
each endpoint p a of a d-edge, draw a line l a through p a that is parallel to the appropriate 
diagonal of C. Let h a be the half-plane defined by l a and not containing the d-edge, and 
set R a = h a f] C. Finally, set R = (n{R a }) n (f){hi}). It is clear that R is convex and that 
any external edge with endpoint in C does not intersect the interior of R. Moreover, by 
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the way we chose the endpoints of the d-edges and s-edges, it is clear that Pq is contained 
in the interior of R. 

3.2.2 Part 2: Putting everything together 

In this part we form the external edges of M = . For each external edge of M' x connecting 
between cells C\ and C2, let p (resp., q) be the point that was chosen as the endpoint in 
Ci (resp., in C2) of the corresponding edge of M = , and match p with q. 

We have already shown that if e\ and e2 are two edges of M=, for which there exists a 
grid cell containing an endpoint of both e\ and e2, then e\ and e2 do not cross each other. 
It remains to verify that a d-edge e of M = connecting between C\ and C2 cannot cause any 
trouble in the cell C through which it passes. Notice that e n C is contained in the danger 
zone in C defined by e (or, more precisely, by the d-edge of M' x corresponding to e). But, 



by Lemma 3.2 (Property 4), this danger zone is empty of points of P, thus e cannot cross 



an internal edge contained in C. Moreover, Lemma 3.2 (Property 2) guarantees that there 



is no d-edge in M= crossing e, and Lemma 3.2 (Property 6) guarantees that there is no 



s-edge crossing e. We conclude that e does not cause any trouble in C. 

Finally, observe that bn(M = ) < 2yWS = 2yl0 • 6n(M x ). This is true since the length 
of a d-edge in M = is at most 2\/25 (i.e., the diagonal of S$(a)), the length of an internal 
edge in M= is at most 4<5 (i.e., the diagonal of a single cell), and the length of a s-edge in 
M = is at most 2\/l05 (i.e., the diagonal of a pair of cells sharing a side). 

The following theorem summarizes the main results of this section. 

Theorem 3.4. Let P be a set of 2n points in the plane. Let Mq (resp., Mj$c) be a 
bottleneck matching (resp., a bottleneck non-crossing matching) of P. Then, 

< 2V10. 

2. One can compute in 0(n 1-5 log ' 5 n) time a non-crossing matching M of P, such that 
bn(M) < 2V^0-bn(M NC )- 
Proof. 1. By applying the algorithm of this section to Mq, we obtain a non-crossing 
matching M, such that 6n(M NC ) < bn(M) < 2 V / 10 • bn(M c ). 

2. Compute Mq in 0(n L5 log ' 5 n) time, using the algorithm of Chang et al. [7j. Then, 
apply the algorithm of this section to Mq to obtain a non-crossing matching M, such that 
bn(M) < 2\/T0 • bn(Mc) < 2\/l0 ■ bn(M^jc). It is easy to see that the time complexity of 
the latter stage is only O(nlogn). ■ 

Remarks. 1. We can improve the approximation ratio to (1 + \^2)^/5 by reducing the 
cell size to 1 + y/2, and picking the endpoints of the d-edges more carefully. 
2. There exists a set P of 2n points in the plane, for which 6ra(MNc) > -^p ■ bn(Mc); see 



Figure 11 



3. It is interesting to note that in the bipartite version, the ratio bn(M^c) /bn(Mc) can 
be linear in n, even if the red and blue points are separated by a line; see Figures [i~2j 
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Figure 11: The length of the solid edges is 1 and the length of the dashed edges is y/85/8. 
The bottleneck (crossing) matching consists of the 8 solid edges, while any bottleneck 
non-crossing matching must use at least one of the dashed edges. 



(a) (b) 

Figure 12: (a) Crossing and (b) non-crossing matching of n blue and n red points separated 
by a line. 

4 Solving Special Cases 

In this section, we consider two special cases of the BNCM problem: when the points of 
P are in convex position, i.e., the points in P form the vertices of a convex polygon, and 
when the points are located on a circle. 

4.1 Matching Points in Convex Position 

Let {pi,p2, ■ ■ ■ ,P2n} denote the vertices of the convex polygon, that is obtained by con- 
necting the points in P, ordered in clockwise-order with an arbitrary first point p\\ see 



Figure 13 Notice that, a non-crossing perfect matching in {pi,j>2, • • • ,P2n} always exists. 
Let M* be an optimal matching of P, i.e., a non-crossing matching with minimum bottle- 
neck. We first observe that, for each edge (pi,Pj) in M*, i + j is odd. According to this 
observation, we define the following weight function 

: if i + j is odd; 
: otherwise. 

Notice also that, for any 1 < i < j < 2n, since the points in {pi,P2, ■ ■ ■ ,P2n} are in 
convex position, so are the points in {pi,Pi+i, ■ ■ ■ ,Pj}- Let M[i,j] denote the bottleneck 
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Figure 13: The convex polygon that is obtained from P, for n = 4. p\ can be only matched 
to the points P2,P4,Pe or p$. 



of an optimal matching between the points in {pi,pi + \, . . . ,pj}. Since M* is a perfect 
matching, each point, particularly p±, is matched in M* . Let pk (for an even k) be the 
point that is matched to p\ in M* . Hence, the bottleneck of M* (i.e., M[l,2n]) is equal 
to max{wi t k, M[2, k — 1], M[/c + 1, 2n]}. Thus, in order to compute M[l, 2n], we compute 
maxjioi fc, M[2, k — l],M[k + l,2n]} for each even k between 2 and 2n, and we take the 
minimum over these values. In general, for every 1 < i < j < 2n, we have 



M[i,j] = 



mm 

=i+l,i+3,. 



max{w it k,M[k + 
max{wi k, M[i + 1, k — 1]} 
max{wi, k ,M[k + l,j],M[i + l,k-l]} 



if k = i + 1 = 
if fc = i + l; 
if k = j; 
otherwise. 



j; 



We can compute M[l, 2n] using dynamic programming. The dynamic programming 
table M has n rows and n columns. Each cell M[i,j] corresponds to a solution of the 
problem for the set {pi,pi + i, . . . ,Pj}, and it can be computed by at most n lookups in the 
table. We fill the table iteratively in such a way insuring that, for each cell M[i,j], all 
the values needed to compute M[i,j] are already computed. To do that, we first fill the 
diagonal M[l, 2], M[2, 3], . . . , M[n — 1, n], and then we fill the above diagonals one by one. 
This leads to solve the problem in 0(n 3 ) time and in 0(n 2 ) space. 

The following theorem summarizes this result. 

Theorem 4.1. Given a set P of2n points in convex position, one can compute a bottleneck 
non-crossing matching in time (9(n 3 ) and in space 0(n 2 ). 



Remark. The same result holds for computing a bottleneck non-crossing matching of 
two (red and blue) sets of points in convex position, and of a set of points on the boundary 
of a simple polygon, with the constraint that edges must not leave the polygon. 
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4.2 Matching Points on a Circle 

Let P = {pi,P2, ■ ■ ■ ,P2n} be a set of points located (in clockwise-order) on a boundary of 
a circle. In this section, we show how to find a bottleneck non-crossing matching of P in 
time 0{n). 

Let M* be an optimal matching of P, i.e., a non-crossing matching with minimum 
bottleneck. We first claim that each edge in M* connects between two consecutive points 
from P, i.e., if (pi,Pj) in M* then j — i = ±1 (assume that p\ = pm+i and p2 n = Po)- To 
see that, assume that there is an edge (pi,Pj) in M* that violates the claim, and consider 
the line that passes through pi and pj. It divides the circle into two arcs pipj and pjpi- 
Assume without loss of generality that pip~j is shorter than pyp%. Then, it is clear that the 
distance between any two points on the arc pip~j is at most PiPj, and, since there is an 
even number of points on pipj, we can match the points pi,pi+i, ■ ■ ■ ,pj-i,pj consecutively 
without affecting the bottleneck of M*. 

Therefore, since there are only two ways to match the points consecutively, we can 
compute an optimal matching of P in linear-time. 

Theorem 4.2. Given a set P of 2n points on a circle, one can compute a bottleneck 
non-crossing matching in time 0(n). 
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